The efficient generation of unstructured control 
volumes in 2D and 3D 



Jacek Leszczynski, Sebastian Pluta 

Technical University of Czestochowa, 
Institute of Mathematics and Computer Science, 
ul. Dabrowskiego 73, 42-200 Czestochowa, Poland 
jale@k2.pcz . czest .pi, pluta@matinf .pcz . czest .pi 



Abstract. Many problems in engineering, chemistry and physics require 
the representation of solutions in complex geometries. In the paper we 
deal with a problem of unstructured mesh generation for the control 
volume method. We propose an algorithm which bases on the spheres 
generation in central points of the control volumes. 



1 Introduction 

An unstructured mesh generator can solve many of the problems associated 
with the structured meshes. The unstructured meshes are suitable for complex 
geometries, especially when we use the control volume method flip} . The control 
volume method on the one hand, is based upon the simple physical principle of a 
flux balance. Nevertheless, the generated meshes have to fulfill some restrictions 
for such method. One of the conditions says than lines connecting the central 
point in the control volume with neighbour central points have to be perpendic- 
ular to the sides of the neighbour volumes. However, the neighbourhood of each 
neighbour or cell in an unstructured mesh must be defined explicitly. This is 
one disadvantage of unstructured meshes because it has to reserve large storage 
of the computer memory. Nevertheless, the advantages like easy for handling 
adaptability in time, ability to generate meshes about arbitrary geometries. We 
can divide such meshes into the cell- vertex control volumes and the cell-centered 
control volumes. In the cell-vertex volumes, all elements containing the relevant 
point are applied as a control volume, so the control volumes overlap. In the cell- 
centered control volumes the elements are subdivided and the control volumes 
do not overlap. 

In this paper we propose a generator for unstructured cell-centered volumes, 
which is useful for the control volume method (lO) . The mesh is constructed for 
a number of points randomly located in a domain and the domain border. 



2 Mathematical background 



We turn our attention in two and three dimensional space. Let us consider a do- 
main Q C R 3 with smooth its boundary dfi. In such domain and boundary we 



generate randomly a set of points Pi(a,b,c), where i — 1..N. The points estab- 
lish cell-centers, around of which we determine convex polygons. The polygons 
are the non-overlapping control volumes which are applied in the control volume 



method 1 1 . We can formulate two groups of basic assumptions necessary for 



the meshes construction as follows: 

1. global conditions: 

- such domain has to be consistent, 

- edges of some polygon have the same lengths to the corresponding edges 
of neighbour polygons, 

2. local conditions: 

- inside a polygon one may found only one random point called the cell- 
center or the point, 

- an edge of the polygon has to be perpendicular to the line connecting two 
points. 

If we fulfill the global assumptions, we may discretize the domain Q through its 
division into the triangles in R 2 space or into the tetrahedrons in R 3 space. Such 
discretization is a standard solution, which one may find in convex geometry 
[Q through Delaunay triangulation |p],p|JT]J|. Regarding to the global and local 
conditions we take into consideration a fact, that we cannot generate polygons 
in R 2 which overlap neighbour polygons. The fact also exist for the generated 
tetrahedrons in R 3 space. 

In this paper we propose a novel algorithm for generation of control volumes. 
We start the algorithm in two dimensions because one can see how it works. We 
consider circles with unknown radiuses which centers are located in the points, 
which are randomly located in the domain fi. Let us assume, that within the 
triangle created by the points the circles only in the one point crosses. The point 
is a vertice of control volumes generated around the point. Following that we 
need to solve the system of equations 

(x-a l ) 2 + (y-b l ) 2 =rl(l = 1,2,3). (1) 

Assuming, that the radiuses r; of the circles are unknown, we obtain a point 
Q(x{r\, T2, 7-3), y(ri, ^2, ^3)) which coordinates are solution of the system ([!]) 



w s -(b2-b 3 )-r 2 1 +(b 1 -b 3 )-r 2 2 -(b 1 -b 2 )-r 2 3 + (a 2 1 +b 2 1 )-(b 2 ~b 3 ) 

X[n, r 2 , r 3 ) - 2-[(oi-aa)-(6i-6 s )-(oi-03)-(6i-63)] + 

(^+^)-(h 1 -& 3 ) + (^+^)-(bl-h2) 

2-[(ai-a 2 )-(6i— 63)— (ai— a3)-(6i— 62)] 
ZT(~ „ „ \ _ (a2-a 3 )-r 2 -(a 1 -a 3 )-r 2 + (a 1 -a 2 )-r 2 -(a 2 +b 2 )-(a 2 -a 3 ) , 

y(ri, r 2 , r 3 ) - 2-l(a 1 -a 2 )-(b 1 -b 3 )-(a 1 -a 3 )-(b 1 -b 2 )\ + 

{ai,+b%) ■(ai-a 3 )-(a§+&i)-(ni-n2) 
~ r 2-[(a 1 -a 2 )-{b 1 -b 3 )-(a 1 -a 3 )-(b 1 -b 2 )] 



(2) 



The indexes in formulas (0) and (||) correspond to the local case. 

In the global case we have random points Pi(a,b, c),i — 1..N which are 
generated inside the domain Q and its boundary dfl. We use Delaunay triangu- 
lation H^pl] that to find neighbourhood of point Pi defined by several points 
PjU,x) ■ ■ ■ Pj(i,k) ■ ■ ■ Pj(i,M(i)) ■ The index k) is a function which establishes 



a point number having neighbourhood with the i-th point. The temporal index 
k varies from 1 to M(i). The M(i) is a number of points corresponding to the 
i-th point. Three points Pi, Pj(i,k) > Pj(i,k+i) define a triangle. When the tempo- 
ral index k + 1 exceeds M(i) then we have PjU^+i) = Pj(i,i)- Moreover, we 
take into consideration a condition j(i,k) ^ i. We extend a definition of for- 
mula (|) for the global case putting (a,, hi), (a j{lM) ,b j{lM) ), (a^ iik+l ), 6 J -( iife+1 )) 
and (fj, Tjt^fAi r j(i,k+i)) respectively into variables with established indexes. In 
our approach, we are looking for the radiuses r^. Taking into consideration of 
such fact we formulate a range of radiuses variation as 

max [(Zjj/j m — r max j(j,fe))] < fj < r max i, (3) 

k 

where L^jak) is a distance between two neighbour points defined as 

= \J («i - a j(i)) 2 + (h ~ b j(i)) 2 - ( 4 ) 
Maximal radius of a circle is formed as 

r m axi = mm [/ii,j(i,fc)] j ( 5 ) 
where is a distance defined in acute triangles as 

^ _ - b j(i,fc+i)) a ~( a j(',fc)^ a j(^fc+i)) b 'l | 

IJ< -*' fe ' V(bj(i,fc)-6j(i,fc + l)) 2 + (aj(i,fc)-ij(i,'« + i)) 2 /g\ 

I l( a j(i,fc)- a j(i.fc + l)) b 3(i,fc + l)-( b j(i,fc)- b 3(i,fc + l)) a 3(i,fc + l) ' ^ ' 

V( b j(i^) _6 3(i,'! + l)) 2 + ( a i(i,'!) _a j<i,'= + l)) 2 

In acute angles of right and obtuse triangles we have 

hij(i,k) — m ^ n [Lij(i^k)i Liju^k+i)] ■ (7) 
k 

The symbol r max j(i ^ represents a maximal radius of the j'-th point being 
a neighbour of the i-th point. Using extended formula (|^) and additional ex- 
pressions (|))-^(@) one can construct an unstructured mesh. Fig. |l| shows an idea 
of the control volume generation using the expression (||) and the limit (||). We 
can observe that the local condition - an edge of the polygon has to be perpen- 
dicular to the line connecting two points - is fulfilled. The general algorithm uses 
some optimisation that to find the points (||) being the vertices of the polygons 
(control volumes). In this paper we neglect the construction of an objective func- 
tion because it is not the subject of our consideration. We applied Rosenbrock's 
method [6] for the local solution. For the global solution we use the soft selection 
method [3| which is a global optimisation method. Prediction of the unstructured 
mesh under such way is not only the one way of the mesh construction. It is pos- 
sible to construct another type of mesh without restriction that expression (||) 
is a solution of eqn. ([!]). But we have to restrict the expression because it 
guaranties right mesh construction for the control volume method JTo| . Let us 
assume that a radius of the circle overlaps other radiuses of the neighbour circles. 



- lines predicting polygenes (control volumes) 

- lines connecting nodes 

- minimal radius of a circle 

- maximal radius of a circle 

- successful area 



Fig. 1. Idea of generation of control volumes. 

Of course we have to restrict the condition (0) . Corresponding to previous and 
present mesh construction we have 

L ij(i,k) - (n + r j(i,k)) < o. (8) 

Following eqn. (||) the point coordinates given by expression (0) are not the 
solution of the system ([!]) but they become the coordinates of the point predicted 
by lines connecting two points of crossing circles. From the other hand, we can 
assume that a radius of the circle does not overlap other radiuses of the neighbour 
circles. In such situation we have 

L i,j(i,k) ~ (n + r j(i,k)) > 0. (9) 

One can say, that conditions (0) and (^|) allow us to obtain different forms of 
unstructured meshes. However, in local solution is possible some combination of 
expressions (||) and @ respectively. 

Similar to previous results we consider a set of random points Pi(a,b,c) in 
i? 3 space. Analogously to the previous considerations, we use spatial Delaunay 
triangulation that to find neighbourhood of point Pj which is defined by several 
points Pj(i.k.i)- The index j(i,k,l) is a function which establishes a point num- 
ber having neighbourhood with the i-th point. The temporal index k describes 



a number of tetrahedrons having the common i-th point. It varies from 1 to 
N(i). The next index I varies from 1 to 3 and it is a point number correspond- 
ing to the fc-th tetrahedron. The four points Pi, Pj(i,k,i) > Pj(i,k,2) > Pj(i,k,3) define 
a fc-th tetrahedron. We search radiuses of spheres which centers are located in 
the points. In each tetrahedron crosses the four radiuses and they create a point 
which is the vertice of the control volume. Solving the following system 



(x - ai f + {y-hf + {z- q) 2 - rf, (I = 1..4), 



(10) 



we can find the local vertices coordinates 



W, 



w z 



x(n, r 2 , r 3 , r 4 ) = — -, y(r u r 2 , r 3 , r 4 ) = — -, z(n, r 2 , r 3 , r 4 ) = — y , (11) 



W 



w 



w 



where 



W 



2(ai - a 2 ) 2{b x - b 2 ) 2(ci - c 2 ) 
2(ai - o 3 ) 2(6i - 6 3 ) 2(ci - c 3 ) 
2(oi - o 4 ) 2(6i - 64) 2(ci - c 4 ) 



Ai 2(6i - 6 2 ) 2(ci - c 2 ) 
A 2 2(6i - 63) 2(ci - cs) 
A 3 2(6i - 64) 2(ci - c 4 ) 



and 



2(oi - 02) Ai 2(ci - c 2 ) 
2(oi - a 3 ) A 2 2(ci - c 3 ) 
2(ai - a 4 ) A 3 2(ci - c 4 ) 



Ai 
A 2 



2(ai - a 2 ) 2(6i 
2(oi - 03) 2(6i 
2(ai - a 4 ) 2(6i 



62) Ax 

63) A 2 
hi) A 3 



A 3 = rj - r 



b\-b\- 
b\-b\- 
-b\-b\ 



-2' 



-3> 



Similar to previous mesh construction, we define lower and upper limits of the 
radiuses variation 



max 

h 



« [{K,j(i,k,l) - r m axj(i,k)j\ <n <r n 



(12) 



The symbol hiju^n represents height of a triangle being the wall of the tetra- 
hedron and for the acute triangles is defined by the following formula 



, _ / (BB 12 - 

n i,j(i,k,i) - y 



CC 3 2-BB 32 -CC 1 2y(AA 32 -CCi2-AA 12 -CC32r+ , 
AAl 2 + BBl 2 +CCf 2 + 



Ram 



AA*+BB*+CC* 



(13) 



where 



AA12 = ^ — ttj(i,fe,j), 
BB12 = h — bj(i t k,i), 
CC12 — Q — Cj( it k,i), 



AA 32 
BB 32 



1 j(i,kJ + l) ~ a j(i,k,l) 1 



Oj(i,k,l+l) 



CC 32 — Cj(i t k,i+i) - c 



°j(i,k,l), 
j(i,k,l) ■ 



In acute angles of right and obtuse spatial triangles we have 

hij(i,k,l) = m i n [Li,j(i,k,l)i ^i,j{i,k,l+l)\ > (14) 

where L^ju^s is a spatial distance defined between two points as 

Li,j(i,k,l) = \J (ai - a j(i,k,i)) 2 + (bi - bj( it k,i)) 2 + (ci ~ Cj(i,fej)) 2 . (15) 

When the temporal index / + 1 exceeds 3 then we put PjU t k,i+i) = Pj(i,k,i)- 
Maximal radius of a sphere is defined as 

r ma xi = nfin [H itj ( iik )] . (16) 

However, the function Hi t jU]A is height of a tetrahedron created by the vertices 
Pi, Pj(i,k,i), Pj(t,k,2), Pj(i,k,3) and for the acute tetrahedron we have 

\A-ai + B-bi + C-Ci + D\ 
H idm + B . 2 + C2 (17) 

where 

A = m 21 • "31 - m 31 ■ n 2 i, B = d 3 i ■ n 2 i - d 21 ■ n 3 %, 

C = d 2 i ■ m 3 i — d 3 i ■ m2i, D = —A ■ aj^^t) - B ■ bj( i>k> i) - C ■ Cj^^t) 

and 

<fei = aj(i,k,2) ~ a j(i,k,i)i m 2i = bj(i t k,2) — bj(i t k,i), n 21 — c^ i k 2 ^ — Cj( i>ktl j, 

d-31 = Oj(i,k,3) ~ a j(i,k,l)i m 31 — bj(i t k,3) — kj(i,fc,l)i "31 = c j(i,k,3) ~ c j(i,k,l)- 

But for the right and obtuse tetrahedrons we have 

Hi,j(i,k) = m } n [Kj(i,k,l)] ■ (18) 

The symbol r max jujA represents a maximal radius of j-th point being a neigh- 
bour of the z-th point. 



3 Example of calculations 

In practical mesh construction when expression (|^) is fulfilled we obtain a unique 
solution. Fig. || shows an example of mesh created under restriction that for- 
mula (||) is a solution of the system ([TJ). Fig. ||a presents triangulation of some 
points in which generated circles. Within the triangle crosses the circles only 
in one point. According to the previous description we draw control volumes 
presented by Fig. ^b. Nevertheless, we can construct a mesh without restriction 
that formula ^ is a solution of the system (|l|). We can apply conditions (||) 
or (||) respectively. In this case we generate a mesh much more easier than in 
previous case. Expressions (||) and (0) extend our considerations for different 
generation of meshes. 




Fig. 2. Example of control volumes taking into consideration that formula (g) 
is a solution of the system (Q): a) random points after triangulation; b) forms of 
control volumes. 



4 Concluding remarks 



In this paper we elaborate a new method of generation of unstructured meshes. 
The algorithm is prepared for the control volume method in two- and three 
dimensional space. For random location of points we use Delaunay triangulation 
before the mesh generation. Our method bases in 3D on circles generation in 
the points. Within one triangle created by three points crosses the circles in one 
point. The point is a vertice of a control volume. We extend our considerations for 
the triangles in which are overlapping or non-overlapping circles. The extension 
allow us to generate meshes much more easier than in previous case. In opposite 
to the Thyssen polygons we can generate control volumes for the right and 
obtuse triangles. However, our method is suitable and easy to use in generation 
of control volumes in three dimensional space. 
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